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^ ; Abstract 

The dynamics of a symmetry breaking phase transition is studied in a ra- 
diation and matter dominated spatially flat FRW cosmology in the large N 
limit of a scalar field theory. The quantum density matrix is evolved from an 
initial state of quasiparticles in thermal equilibrium at a temperature higher 
than the critical. The cosmological expansion decreases the temperature and 
triggers the phase transition. We identify three different time scales: an 
early regime dominated by linear instabilities and the exponential growth of 
long-wavelength fluctuations, an intermediate scale when the field fluctuations 
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Q^' probe the broken symmetry states and an asymptotic scale wherein a scaling 

O ■ regime emerges for modes of wavelength comparable to or larger than the 

horizon. The scaling regime is characterized by a dynamical physical correla- 
tion length S,ph ys = dn{t) with dn(t) the size of the causal horizon, thus there 
is one correlated region per causal horizon. Inside these correlated regions 
the field fluctuations sample the broken symmetry states. The amplitude of 
the long-wavelength fluctuations becomes non-perturbatively large due to the 
early times instabilities and a semiclassical but stochastic description emerges 
in the asymptotic regime. In the scaling regime, the power spectrum is peaked 
at zero momentum revealing the onset of a Bose-Einstein condensate. The 
scaling solution results in that the equation of state of the scalar fields is the 
same as that of the background fluid. This implies a Harrison-Zeldovich spec- 
trum of scalar density perturbations for long- wavelengths. We discuss the 
corrections to scaling as well as the universality of the scaling solution and 
the differences and similarities with the classical non-linear sigma model. 
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I. INTRODUCTION, MOTIVATION AND SUMMARY 



Symmetry breaking phase transitions play a fundamental role in the Standard Model 
of particle physics and its generalizations and are conjectured to have ocurred in the early 
Universe. Some versions of inflationary models invoke a symmetry breaking phase transition 
to generate the inflationary stage [0J and the process of field ordering after the phase transi- 
tion has been argued to produce density fluctuations that could seed galaxy formation [jlJQ . 
Certainly, symmetry breaking phase transitions in the early Universe are very important at 
the electroweak scale and at much lower temperatures for the hadronization and chiral phase 
transitions in QCD. At a more speculative level, late time phase transitions in the matter 
dominated era had been proposed || with the possibility that pseudo Nambu-Goldstone 
bosons may contribute to the dark matter component and influence large scale structure 
formation. 

Typically, symmetry breaking phase transitions are studied by means of an effective po- 
tential which is the free energy density for a space-time constant expectation value of the 
order parameter. Although this formulation provides a good qualitative understanding of 
the equilibrium properties it is unable to address the important issue of the dynamics of 
symmetry breaking and the non-equilibrium process of phase ordering and phase separa- 
tion. The simple picture of the order parameter rolling down the potential hill although 
intuitively appealing does not capture the correct dynamics. A consistent non-equilibrium 
field theoretical treatment reveals that the quantum fluctuations are extremely important 
and enhanced by the instabilities that are the hallmark of the phase transition Q. 

The picture that emerges from consistent studies of the dynamics of phase transitions in 
Minkowski space time 0-0] is remarkably similar to phase ordering kinetics in a symmetry 
breaking phase transition in condensed matter systems |^-[T0[]. 

Recent results in quantum field theory |],[?[] in Minkowski space-time reveal the following 
picture: when the system is cooled below the critical temperature on time scales shorter 
than the relaxation time of long-wavelength fluctuations, the process of phase separation 
begins via linear spinodal instabilities and the exponential growth of these long-wavelength 
fluctuations ||^,[7|. A new dynamical time scale emerges at which the back-reaction of 
these fluctuations becomes of the order of the tree level term and the linear instabilities 
shut-off. At this time scale, the fluctuations of the field sample the broken symmetry equi- 
librium states (minima of the free energy) which in turn implies that the amplitude of the 
long-wavelength fluctuations becomes non-perturbatively large f4j||0] and consistent non- 
perturbative methods are necessary to study the real-time dynamics of the phase transition. 
When the backreaction from the field fluctuations becomes comparable to the tree level 
terms in the equations of motion the non-linearities become very important. In this late 
time regime a dynamical length scale emerges £(£) = t which is a consequence of causality 
and determines the size of the correlated regions, the emergence of a scaling behavior and 
the onset of a novel non-equilibrium condensate |6|-p[pT|. Recently a very thorough study 
of the non-equilibrium time evolution has been provided in finite volume (finite momentum 
resolution) with non-perturbative techniques providing a consistent picture of the dynamics 
of phase separation ][iT| . 

In cosmology the non-equilibrium dynamics of a symmetry breaking phase transition 
has been studied within the context of inflationary scenarios including quantum backreac- 
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tion from the matter fields in the metric dynamics, and non-equilibrium aspects of particle 
production and reheating including quantum backreaction in the large N limit had been 
studied in Friedmann-Robertson- Walker (FRW) cosmologies in a fixed background fll"3| , |I4 



We investigate here the dynamics of a phase transition and phase ordering kinetics in a 
radiation and in a matter dominated FRW background metric from a quantum field the- 
ory formulation that includes self-consistently the quantum and thermal fluctuations in the 
dynamics. 

Phase ordering dynamics in a FRW background has been previously studied in the large 
N limit of the classical non-linear sigma model in []I5,IE|. Detailed numerical and analytical 



studies of the classical non-linear sigma model revealed the emergence of a scaling solution 
of the classical field equations that reflect the dynamics of Goldstone bosons in the broken 
symmetry state |TI|. The scaling solution allowed to study the cosmological perturbations 
seeded by the fluctuations of the Goldstone scaling field fl7||T6ll and provides a solvable 
example to study the evolution of causal cosmological perturbations. The importance of a 
scaling solution lies in the fact that for causal perturbations the dynamical range of time 
scales between the phase transition and matter-radiation decoupling is very large to be 
tractable numerically, the scaling property of correlation function thus allows to extrapolate 
them to arbitrarily long time scales P,CT- 

The non-linear sigma model is a theory that describes the dynamics of Goldstone bosons 
in a broken symmetry state and as such the field is constrained to the vacuum manifold. 
Hence this model does not allow to study the question of the dynamics of symmetry breaking 
phase transitions and the non-equilibrium aspects of phase separation and ordering beginning 
from an initial high temperature state in the unbroken symmetry phase and cooling through 
the phase transition until correlated regions of broken symmetry are formed and grow. 

The goals: The goals and focus of this article are to study the important issue of the 
dynamics of the phase transition which simply cannot be addressed within the context of the 
non-linear sigma model. Namely, i) the non-equilibrium dynamics of a symmetry breaking 
phase transition from an initial high temperature unbroken phase to a final low temperature 
broken symmetry phase in a radiation or matter dominated background FRW cosmology, 
ii) to follow the process of phase separation and ordering directly in real time from the 
early time instabilities to the formation of correlated regions, iii) to study the emergence 
of a scaling regime as a dynamical consequence of the phase transition and iv) implications 
for the equation of state of the matter described by the quantum scalar field. Furthermore 
we analyze the deviations from scaling and discuss the necessity for scaling violations for 
self-consistency. The description of cosmological perturbations is outside the scope of this 
paper. We keep our focus on the description of the non-equilibrium dynamics of the phase 
transition and the resulting consequences. 

The model: As argued above a consistent study of the dynamics of a phase transition 
and the process of phase separation requires a non-perturbative framework. There are 
very few non-perturbative frameworks that are i) renormalizable ii) maintain all of the 
conservation laws, iii) lend themselves to a detailed analytical or numerical study and iv) 
can be consistently improved. We thus study the quantum linear sigma model of a scalar 
field in the vector representation of O(N) in the leading order in the large N limit in a fixed 
spatially flat radiation and or matter dominated FRW cosmological background. There are 
several noteworthy features of this model that deserve comparison with the classical non- 
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linear sigma model studied in detail in ||T5|JT6[| : i) we consider the full quantum field theory as 
described by a time dependent quantum density matrix prepared initially in equilibrium at 
high temperature, ii) unlike the non-linear sigma model, this is a renormalizable quantum 
field theory and the ultraviolet divergences are absorbed in renormalizations of masses, 
couplings and field amplitudes, iii) the linear sigma model allows to go from the unbroken 
to the broken symmetry phase, the non-linear model is only defined in the broken symmetry 
phase. 

The strategy: The main ingredient in the description of the non-equilibrium dynamics 
is the time evolution of an initially prepared quantum density matrix. We begin our analysis 
by implementing the large N limit, which to leading order leads to a self-consistent quadratic 
Hamiltonian with effective frequencies that depend on time through the background scale 
factor as well as through the self-consistent mean field. The density matrix is taken to be 
Gaussian consistent with the leading order in the large N and its initial form describes a 
state in which the self-consistent quasiparticles are in local thermal equilibrium at an initial 
temperature Tj > T c with T c the critical temperature. We will not justify the choice of an 
initial thermal state for the self-consistent quasiparticles, and we will simply assume that 
such initial state is physically reasonable and is an acceptable description of the system prior 
to the phase transition. This initial density matrix evolves in time through the Liouville 
equation which has an exact solution in the leading order in the large N limit. The advan- 
tage of our formulation in terms of the quantum density matrix is that it leads to a clear 
description of the emergence of (semi) classical but stochastic description. The expansion of 
the cosmological background results in adiabatic cooling that triggers the phase transition 
resulting in long-wavelength instabilities which are the hallmark of the process of phase sep- 
aration and ordering and completely dominate the early time evolution of correlations. The 
early time dynamics can be understood analytically since for weak coupling the backreac- 
tion is perturbatively small and the equations are essentially linear. The intermediate and 
asymptotic dynamics are studied analytically and numerically both in radiation and matter 
dominated FRW backgrounds. 

Summary of results: The main results of this article can be summarized as follows: 

• The dynamics of the phase transition reveals several different time scales: an early 
time scale dominated by the exponential growth of long-wavelength fluctuations as a 
consequence of the instabilities associated with the phase transition. An intermediate 
time scale at which the self-consistent backreaction from the quantum and thermal 
fluctuations begins to compete with the tree-level terms in the equations of motion 
and signals the onset of non-perturbative and non-linear dynamics. An asymptotic 
time scale and the emergence of a scaling solution for fluctuations with wavelengths of 
the order of or larger than the causal horizon. Corrections to scaling are an unavoidable 
consequence of the dynamics and they are important for self-consistency, their form is 
uniquely determined by the background. 

• In the scaling regime there emerges a dynamical length scale £ p hys{t) — dnif) with 
du{t) the size of the causal horizon. This length scale determines the size of the cor- 
related regions, inside which the field fluctuations probe the broken symmetry states. 
There is one such correlated region per causal horizon. This is a microscopic, quantum 
field theoretical justification of Kibble's original proposal PJI7||. The power spectrum 
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is peaked at superhorizon wavelengths and signals the onset of a non-equilibrium con- 
densate. 

• After the intermediate time regime the amplitude of long-wavelength fluctuations 
become non-perturbatively large and the phases freeze out. A semiclassical but 
stochastic description emerges and field configurations with large amplitudes and long- 
wavelengths are represented in the ensemble (density matrix) with non-negligible prob- 
ability. 

• As a consequence of the scaling solution, the equation of state for the fluid of the field 
fluctuations is the same as that of the background. As a consequence, the density 
fluctuations 5p/ Pbackground are time independent. Therefore, superhorizon models will 
exhibit a Harrison-Zeldovich spectrum. 

• The universality of the scaling solution and the corrections to scaling are discussed in 
detail. 

The article is organized as follows: In section II we set up the model, implement the 
large N limit to leading order and solve the equations that determine the time evolution of 
the initially prepared density matrix. In section III we study analytically and numerically 
the dynamics of symmetry breaking, establishing the different time scales, the emergence 
of scaling and of a semiclassical description. In section IV the consequences of the scaling 
solution are analyzed: the dynamical length scale and the correlated regions and the equation 
of state of the scalar fields fluctuations. Section V provides a comparison between our results 
and those obtained for the classical O(N) non-linear sigma model. Section VI summarizes 
our conclusions. 



II. DENSITY MATRIX EVOLUTION IN THE LARGE N LIMIT 

We start by setting up the formulation of the time evolution of the quantum density 
matrix in a spatially flat FRW cosmology with a fixed background metric. In comoving 
coordinates the metric element is given by 

ds 2 = dt 2 -a 2 (t) dx 2 (2.1) 

with 




and p(t) is a fixed background energy density which is taken to be a combination of radiation 
and matter with p l R , p l M the energy densities in radiation and matter respectively at the 
initial time t = t with a(t ) = 1. 

We consider a theory of N scalar fields in the vector representation of O(N) in the leading 
order in the large N limit. The action is given by 
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(2.5) 
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with 1Z the Ricci scalar and £o has been introduced with the purpose of renormalization. 
Since the time evolution of the quantum density matrix is determined by Liouville's equation, 
we need the Hamiltonian which requires the canonical momentum conjugate to <$>(x,t) 



The Hamiltonian becomes 



H(t) 



(Fx 



U(x,t) = a 3 (t) $(x,t) 



n 2 (f,t) , a W ( v$)2 + G 3 (t)v(< 5 ) 



2 a 3 (t) 



+ 



(2.7) 



(2.8) 



In the Schrodinger representation (at an arbitrary time t), the canonical momentum is 
represented as 

n a (f) = -i — — ;a = l,---N 

and the functional density matrix p with matrix elements in the Schrodinger representa- 
tion p[$ a (7), $ 6 (7); t] obeys the Liouville equation, which in the Schrodinger representation 
becomes the functional differential equation |12[] 



H 



H 



8®Q 



p[$ a O,<l> b (r);t] 



(2.9) 



To leading order the large N limit can be implemented by the following Hartree-like factor- 
ization 0,0] (for an alternative formulation of the large N limit see [ [18|j ) 



($ • $) 2 -> 2 ($•$)$•$-($• $) 2 ; ($•$)= iV($ a $ a ) (no sum over a) 



(2.10) 



where the expectation value is in the time evolved density matrix p(t) which is the solution 
of the Liouville equation above. 

It is convenient to introduce the spatial Fourier transform of the fields as 



$(£, t) 



Ak-x 



(2.11) 



with Q the spatial volume, and a similar expansion for the canonical momentum U(x,t). 
The Hamiltonian becomes 



k 



a 3 (t) 



2a 3 (t) nrn -* + 2 



W 2 k {t) $ r $ 



w 2 k {t) = -ml + i n + 



k 2 

a?{t) 



d 3 k 
(27p 



($|<I>%) (£) (no sum over a) 



(2.12) 
(2.13) 
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This Hamiltonian describes a set of infinitely many harmonic oscillators, that are only 
coupled through the self-consistent condition in the frequencies (|2.13|) . Since the effective 
Hamiltonian is quadratic in terms of the self-consistent frequencies, we propose the following 
Gaussian Ansatz for the functional density matrix elements in the Schrodinger representation 



p[$, <£, t] = ]J A4(t) exp f-Ml $ fc • $_ fc - ^ j fe ■ l_ fc - B k (t) 4 ■ l_ fc ] (2.14) 
k ^ > 

This form of the density matrix is dictated by the hermiticity condition p^[$, $, t] = 

p*[$, $, t}; as a result of this condition B k (t) is real. The kernel B k (t) determines the 
amount of mixing in the density matrix. If B k = 0, the density matrix corresponds to a 
pure state because it is a direct product of a state vector times its adjoint conjugate. 
The Liouville equation ( pT9| ) becomes fT2" | 




p[$,$,t] (2.15) 



The equations for the kernels in the density matrix (A^, B^) are obtained by comparing 
similar powers of on both sides of the above equation. We obtain the following equations 
for the covariances: 

$=2^*-*> (2 - 16) 
iA k = ^^-a\t)Wt{t) (2.17) 

iB k = ^-{A k -Al) (2.18) 

The equation for B k (t) reflects the fact that an initial pure state with B k (0) = remains 
pure under time evolution. From these equations it becomes clear that the only independent 
equation is that for A k (t). In particular, writing A k (t) = An )k (t) + iAj jk (t) we find that the 
ratios A Rk (t) / B k (t) and Af k (t)/[A Rk (t) + B k (t)\i are constant in time, the latter being a 
consequence of unitary time evolution of the density matrix. Exploiting the proportionality 
between A^ k (t) and B k (t) we introduce a new variable fT2|| 

A k (t)=A Rjk (t)+iA Iik (t) 

by defining 



A R>k {t) = A R , k (t) cothe k (2.19) 

B k (t) = (2.20) 
sinn Bfc 

A Iyk (t) = A Itk (t) (2.21) 
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This new variable obeys the following Ricatti equation of motion 



>^dt) = ^M-a\t)Wl{t) 



a 3 (t) 

This equation can be linearized by defining 



A k (t) = -ia\t) 



and the mode functions (pk{t) obey the following equation of motion 

(p h {t) + 3 ^| Mt) + W 2 k {t) Vk (t) = 
a[t) 



(2.22) 



(2.23) 



(2.24) 



The equations of motion ( 2.24|) are recognized as the Heisenberg equations of motion for $&(£) 
obtained from the Hamiltonian ( [2. 12 ). This observation leads to a rather clear interpretation 
of the mode functions ipkit) as a basis for the expansion of the Heisenberg field operator 
$£(*), i.e. 



V2 



iz<p k {t) - <p* k (t) 



(2.25) 



where the annihilation and creation operators ; a_- are independent of time in the 
Heisenberg picture and define a Fock representation. We will normalize the mode functions 
ipkif) so that the Wronskian is given by 



-2i 



a(to) 
a{t) 



(2.26) 



where to is the initial time at which the density matrix is prepared. Hence we find 



m *«->(t) 



I^M coth e, 

2a 3 (t ) 2 



Va = 1, 



N 



(2.27) 



Without loss of generality, we choose a (to) = 1 which can always be done by a simple 
rescaling of lengths. Furthermore, we will choose 0^ to reflect a thermal density matrix at 
the time t in terms of the self-consistent frequencies which are discussed below. 



A. Conformal time analysis, initial conditions and renormalization 

It is convenient to change variables to conformal time rj defined as 

ft A+i 

V= -7TK 5 V(t = to) = , (2.28) 
Jt a(r) 

which is chosen to vanish at the initial time. 

The scale factor in conformal time C(rj) and conformally rescaled mode functions fk{ff) 
are given by 
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C( V ) = a(t( V )) ■ C(0) 
fk(v) = 



(2.29) 
(2.30) 



The conformally rescaled mode functions obey the Schrodinger-like differential equation 
d 2 



. ]+ k 2 + C 2 ( V ) M 2 ( v ) 



dr\ 

M 2 ( V ) = -m 2 + ^ 
C"{rj) 



1, 
6' 



fk(v) 



o 

^T 2 



fe a/c , u k coth — 
C 2 (?]) 2 



ft(r?) = 6 



C 3 (r?) 



(2-31) 
(2.32) 
(2.33) 



where primes denote derivative with respect to conformal time. We will choose the following 
initial conditions on the mode functions 

1 

(2.35) 



A(o) 



[2.34) 



oo k = ^k 2 + M 2 (Q) 



in this manner the mode functions at the initial time represent positive frequency (particle) 
modes. We now choose 



ft - Uk 



^2.36) 



so that the initial density matrix at t = to] rj = describes a statistical ensemble in lo- 
cal thermal equilibrium with an initial temperature for the conformal modes describing 
(quasi) particles of self-consistent frequencies u k . We will not try to justify this choice of 
initial state in local thermal equilibrium and simply assume that it provides a physically 
reasonable description of the state prior to the phase transition. 

m 



The renormalization aspects had already been studied in detail in references 



with the result that the quadratic and logarithmic divergences in terms of an ultraviolet 
cutoff A can be absorbed in mass, coupling and conformal coupling renormalization Wq — > 
m Ri — *■ ^Ri £o £,r- We refer the reader to those references for details and highlight that 
the main result of the renormalization program is that 



with 



(*P 2 )r 



6{k 



M 2 ( v ) = -m 2 R + ^ R 



Ir + J 

d 3 k \u k | f k (rj) 



A; 



R 



4k 3 



m 



J 



(2tt) 3 [ 2 C 2 (t]) 2 k C 2 (r]) 
K) [ K ( C'(0) Y 

d 3 k | f k (rj) | 2 u k 
(2tt) 3 C 2 (r]) exp^-1 



(2.37) 
(2.38) 

(2.39) 
(2.40) 
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with K an arbitrary renormalization scale [ I2[ . The explicit relation between bare and 
renormalized parameters can be found in [12,13. 19]. The finite temperature contribution 
J can be written by separating in the integral the large k behavior of the mode functions 
fk(v) ~ e~ rkri /^/k for k >> C{rf) M(rj). Assuming that the initial temperature Tj >> M(0) 
we can separate the contribution from the large wave- vectors k > T » M(0) as in the 
hard-thermal loop approximation in thermal field theory |20| , |2l|| 
of J has been obtained in reference |T2] 
thermal loop limit 



The high temperature limit 
the leading contribution is obtained in the hard 



J, 



HTL 



5kM 

12 



T eff(v) = Ti 



C(0) 

c(v) 



Teff(0) 



(2.41) 



This term is completely determined the very short wavelength modes k > T ^> M(0) 
and reflects the adiabatic cooling of the short wavelength modes in the initial state by the 
cosmological expansion. 

Introducing the critical temperature as 



^2.42) 



we find that the hard thermal loop contribution Jhtl combines with the renormalized mass 
term -m\ in M 2 (r]) (|2~32|) to yield 




A 



- m R + '-ttJhtl(v) 



^2 



- 1 



^2.43) 



Clearly the leading order finite temperature contribution, i.e. the hard thermal loop limit 
gives rise to an effective i] dependent mass squared term that changes sign when the effective 
temperature falls below the critical, thus triggering the phase transition. 

There are sub-leading corrections in the high temperature limit that have been studied 



in detail in reference 



The term linear in the initial temperature corresponds to the 



classical contribution for uo k T,, i.e. k Tj which is given by 



d?k | f k ( v ) | 
(2tt) 3 C 2 (r]) 



[2.U) 



The contributions from Jhtl(v) an d Jciiv) are the most important ones in the weak coupling 
limit Afl < 1 as can be seen as follows. The leading order correction Jhtl{v) determines 
the effective squared mass, which is positive for T e ff(r]) > T c and becomes negative for 
T e ff(v) < T c thus triggering the phase transition. Once the effective temperature falls 
below the critical, long-wavelength instabilities begin growing exponentially and the long- 
wavelength modes acquire non-perturbatively large amplitudes as argued in the introduction 
and studied explicitly in detail below. These long- wavelength modes with k << Tj determine 
the contribution J c i{t}), which is the dominant contribution from the infrared sector in the 
high temperature limit. In terms of the critical temperature (|2.42 ) and the ratio Tj/T c > 1 
the classical contribution leads to 



A 



R 



Jcliv) = m R 



y/24A 
4vr 2 



R 



k alk 



fM I 
c 2 (v) 



(2.45) 
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where k is a cutoff that determines the maximum wavevector that will be unstable. Numeri- 
cally it is found that the integral becomes insensitive to the choice of this cutoff for k > tur. 
The next contribution is logarithmic in temperature W2^ and is of order A^ln(A^), therefore 



subleading in the high temperature and weak coupling limit. Furthermore, the zero tem- 
perature part is finite after renormalization and subleading in the weak coupling limit since 
it is of O(Xr) as compared to the (9(1); C(a/Ar) for the 0(Tf); 0(Ti) contributions from 
Jhtl'i Jd respectively. Hence the zero temperature part that remains after renormalization 
will be neglected. 

We finally obtain the equations of motion in leading order in the high temperature and 
weak coupling expansion after the following steps 

• Rescaling all dimensionful quantities, k,rj, fk(rj),Uk, etc. in terms of the only dimen- 
sionful parameter tur. Effectively this amounts to setting itlr = 1 and all dimensionful 
parameters are understood in units of uir. 

• After renormalization of couplings and mass we neglect the zero temperature contri- 
bution, which is higher order in the weak coupling Xr and keep only the leading order 
terms Jhtl(v) an d Jci(v) neglecting higher order contributions which are subleading 
in the high initial temperature and weak coupling limits. 

• We choose the case of minimal coupling £ R = 0. 

In summary, the equations of motion for the mode functions and the initial conditions are 

d 2 



^- 2 +q 2 + C 2 ( V )M 2 ( V ) 



fM = (2.46) 



M 2 (,) = ^-l-gM + ^) (2.47) 
^ = Jo q dq 5 9 = ^^T C (2 - 48) 

W) = - ; f'M = -w* fM (2-49) 



W + ^-l (2.50) 



where we have set the cutoff k = itlr = 1 which will be justified numerically below, and 
neglected perturbatively small corrections of 0(\/Ar) in uj q . The term gY,(rj) describes the 
back-reaction of the quantum and thermal fluctuations. 

With the background energy density determined by the combination of matter and ra- 
diation as given in eq. (|2.3p a simple expression is obtained for C(rj) by introducing the 
quantities 



T2 



Pr ' i + n 3M Pl 

where we have used a(t ) = 1. A straightforward integration leads to 



(2.51) 
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. , H 2 rj o f 1 + H; r\ for Radiation dominated 

C r? = 1 + Hi 7i + ' \ r/ 2 = <^ ^2 „ 2.52 

V /; ' 4(1 + n) ' {l + Hri + ^ri 2 for Matter dominated V ; 

Rescaling Hi and rj in terms of m# we choose the value Hi = 1/2 for convenience. This 
value is motivated by a choice of an initial radiation energy density ~ (10 16 ) 4 (Gev) 4 with a 
mass scale ~ 10 12 Gev. Other values of Hi lead to a proper re-scaling of the time scales and 
therefore a quantitative change in the results, but the qualitative features described below 
are robust. 



III. DYNAMICS OF SYMMETRY BREAKING: 

We begin our study of the dynamics of symmetry breaking by focusing on the case of a 
radiation dominated FRW cosmology, i.e. we set r = in Q2.52j ) and as explained above we 
have chosen Hi = 1/2 for convenience. The dynamics is completely determined by the set 
of equations (|2~46|^30|) with C(rj) = 1 + rj/2. 



Before embarking on a numerical study of these equations we can obtain insight into 
the dynamics by neglecting the backreaction term gT,(rj) in ([2.471) . This approximation is 
valid in the weak coupling limit and for early times when #£(77 ) <C 1. During the time 
scale for which T e jf{rj) > T c the effective squared mass is positive and the mode functions 
are oscillatory functions of conformal time. When the effective temperature falls below 
the critical, the effective squared mass becomes negative and long-wavelength modes with 
momentum q 2 < C 2 (rj) M 2 (rj) become unstable and begin to grow almost exponentially, 
obviously the growth rate is the largest for the smaller q. We can obtain an estimate of the 
behavior of the long-wavelength modes after the temperature falls below the critical via a 
WKB expansion which after approximating M 2 (rj) w — 1 — C"(rj)/C 3 (rj) leads to 

e J v C(r,')\M(rf)\d V ' 
fq~o(ri) / (3- 1 ) 

y/C( V )\M( V )\ 

which reflect an exponential growth of long-wavelength fluctuations as a result of linear 
instabilities of the long-wavelength modes below the critical temperature. This behavior in 
turn entails that the backreaction term also grows exponentially as 

e 2 pC(r,')\M(r,')\d V > 

9m ~ 9 clS (3 ' 2) 

and becomes non-perturbatively large at a time scale r] NL when it begins to compete with 
the tree level term, i.e. when 

gT;(r] NL ) ~ 1 (3.3) 

The condition ( |3.3| ) determines the time scale 77 nl at which the non-linearities become 
important and must be treated non-perturbatively. This condition has a simple physical 
interpretation: it implies that the amplitude of the long wavelength fluctuations probe the 
broken symmetry states, i.e. since 

gZ( V ) « ^(^ a ) R 
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restoring rriR as the scale unit we see that 

with $q being the minimum of the renormalized tree level potential. Once gH ~ ra^ the 
backreaction compensates for the negative squared mass and the long-wavelength instabil- 
ities shut-off. This picture is confirmed by a detailed numerical analysis of the equations 
of motion ( 2.46|) -( p.50 ) with C{rj) = 1 + f]/2 for a radiation dominated FRW cosmology 



(r = 0). 

Fig. H displays the effective squared mass M 2 (rj) given by eq. ( |2.47 ) for a radiation dom- 



inated (r = 0) cosmology with the value of the parameters Tj/T c = 3 ; g — 10~ 5 ; Hi = 1/2 
and all dimensionful quantities had been rescaled by rriR. The effective squared mass be- 
gins with a positive value at rj = since Tj > T c but diminishes as a consequence of the 
adiabatic cooling through the cosmological expansion. The back reaction gH(rj) is displayed 
in Fig. it is initially perturbatively small during the stage in which T e ff(rf) > T c . Once 
the effective temperature falls below the critical, long- wavelength fluctuations begin to grow 
approximately as in eq. ([3.1|) as a result of the linear instabilities and the back reaction term 
begins to grow exponentially and to compete with the tree level term. The rapid expo- 
nential growth of the back-reaction terms results in an overshooting and the effective mass 
M 2 (rj) becomes positive. When this happens the mode functions become again oscillatory 
and their amplitude diminishes causing a damped oscillatory behavior in gT,(rj) and conse- 
quently in M 2 (rj). However, in the Schrodinger-like equations for the mode functions ( |2.46| ) 
the mass term is C 2 (t])M 2 (i]), thus while M 2 (rj) becomes positive and oscillatory and even- 
tually damps out by dephasing of the oscillations the prefactor C 2 (rj) makes this effective 
mass to grow during some period of time. The effective mass C 2 (r])M 2 (r)) for the mode 
functions f q {rj) is displayed in Fig. § which clearly shows the initial cooling, the growth of 
the backreaction and the competition between the damping of the oscillations in M 2 {rf) and 
the growth of the scale factor. This figure also reveals the striking asymptotic behavior that 
C 2 (i])M 2 (r]) — > as rj — > oo. A detailed numerical analysis reveals that asymptotically at 
long time C 2 (r])M 2 (r)) — > 0(l/r] 2 ), a result that will be understood analytically below. The 
combination of figures ([1]-|3|) reveals three different time scales after the phase transition, i.e. 
when T e ff < T c 

• Early time (rj < ^tvl): this stage is dominated by linear instabilities that result in the 
exponential growth of long-wavelength modes. The back-reaction contribution gTi(r]) 
grows exponentially but remains perturbatively small. 

• Intermediate time (r] ~ t]nl): during this stage the back-reaction begins to be com- 
parable to the tree level contribution to the mass, i.e. g'E(^) ~ m\, the back-reaction 
begins to shut-off the linear instabilities but overshoots resulting in an oscillatory 
effective mass. This stage determines the onset of non-linear evolution since the back- 
reaction is of the same order as the tree level term. The damping of M 2 (r) competes 
with the growth of the scale factor C 2 (r)) resulting in that the effective squared mass 
C 2 (r))M 2 (i]) still grows. 

• Asymptotic regime (77 ^> tjnl): This regime is fully non-linear and the detailed nu- 
merical analysis reveals that the effective mass C 2 (r])M 2 (r]) vanishes asymptotically. 
Numerically we find that in this regime 
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CW(J)) «. _(*±^±1) fcr0( , ) = (2)' (3.4) 

with p = 1,2 for radiation and matter dominated respectively. The result ( |3.4|) implies 
a very delicate mechanism of cancellation between the back-reaction and tree level 
terms, which will be understood analytically below. 



A. Asymptotic dynamics and emergence of scaling: 



From the expression for M 2 (rj) given by eq. (|2.47 ) the vanishing of C 2 (r])M 2 (rj) in the 
asymptotic region, when T e ff(rj) <C T c leads to the following sum rule in the asymptotic 
regime rj > r] NL ^ 1 (in units of ttir) 



9 / q 



fM \ 2 = c 2 ( v ) 



2p 



(3.5) 



where p = 1 for radiation dominated and p = 2 for matter dominated. 

The constraint that the integral on the left hand side must lead to a power law suggests 
the following scaling ansatz for the mode functions 



\fM\ 2 = v a \Hx)\ 2 ; x = qv 



(3.6) 



Comparing the powers and the coefficients leads to the following constraints 

a = 2p + 3 (3.7) 



g x 2 dx \F{x)\ 2 = 2- 2p 
Jo 



(3.8) 



Assuming that the scaling function T{x) is regular at x — 0, the scaling ansatz ( |3.6| ) leads 
to a remarkable conclusion: the q = mode function is of the form 



fo(v) = ^(O) 



(3.9) 



leading to the result 



fS{v) _ (2p + 3)(2p+l) 



fo(v) 



-15/477 2 for Radiation dominated 
-35 /Ai] 2 for Matter dominated 



(3.10) 



which is precisely the numerical result ( |3.4|) . The result ( |3.10|) in turn leads to the following 
equations of motion for the q 7^ mode functions 



d 2 



drf Arj z 
with the general solutions in terms of the scaling variable x = qn 



(2p + 3)(2p + i; 



fM = 



(3.11) 
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fqOn) = A a r/2 — h £L —N 2 (x) for Radiation dominated (3-12) 

x 2 r]2 

f q (r}) = A q rfi ^ft- + B g ^jN 3 (x) for Matter dominated (3.13) 

X 6 772 

Where A g , 5 g are numerical constants and J p+ i(x) ; N p+ i(x) are Bessel and Neumann 
functions respectively. The constant Wronskian between f q (rj) and its complex conjugate 
determines that neither A q nor B q could vanish. However in the asymptotic limit 77 — > 00 
and at fixed x, two important simplifications occur: i) the term proportional to N p+ i(x) 
becomes subleading and ii) q = x/r] — > and A q — > A , thus in the asymptotic regime 
for very large 77 and fixed x we find that the asymptotic behavior of the mode functions is 
dominated by the J p+ i(x) contributions with a constant coefficient Aq, i.e. asymptotically 
for rj 3> 1, x fixed, the leading contribution to the solutions are of the scaling form 

f q (n) r? = 1 A jj§ + O \\\ for Radiation dominated (3.14) 

x \ V 2 ) 

f q (r]) = A r]i + O [ 4- J for Matter dominated (3.15) 

x \ V 2j 

In this asymptotic region, the absolute value of the coefficient Aq is completely determined 
by the constraint (|3.8|) resulting from the sum rule (|3.5|) and the integral 



x2p [J P+ 1{X)\ 22p+1 T[2p + | } r(p + 1 )2 



We find, 



2 J i|| = 2.94524. . . for Radiation dominated 
= I W# = !7-39534 .../g for Matter dominated ^ 3 ' 16 ^ 

Figs. fU?] display p r/~ 5 | f q (r] = 240) | 2 ; gr/~ 5 |/ 9 (r/ = 400) | 2 and g\A J 2 (x)/x 2 \ 2 vs. a; respec- 
tively for a radiation dominated cosmology with g = 10 -5 ; Tj/T c = 3. These figures are 
indistinguishable from each other, furthermore we have checked numerically that 77" 2 f q=0 (j]) 
approaches a constant asymptotically both for the real and the imaginary part in a radiation 
dominated FRW cosmology. Similar results had been obtained numerically for the case of 
matter domination, in particular confirming the scaling behavior ( |3.15| ) with the appropriate 
coefficient given by ( |3.16| ). In this case also both the real and imaginary part of rj~2f q=0 (rj) 
approach a constant asymptotically and the effective mass in both cases oscillates with small 
amplitude around a mean value given by ( p.lOj ) as is shown explicitly in fig. (H). 

Thus a detailed numerical integration of the equations of motion confirms the scaling 
ansatz and leads to the following conclusions: 

• Asymptotically the effective mass C 2 (77) M 2 (if) n = 1 —Qp+^2s±H_ The vanishing of the 
effective mass leads to the sum rule ( p.5| ) which suggests the scaling ansatz ( |3.(p with 
the power law given by ( |3.7|) ). The self-consistent solution in the asymptotic regime is 
given by (|3.12|) for radiation dominated and (|3.13|) for matter dominated. Neither of 
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the coefficients A q ; B q vanishes as a consequence of the constancy of the Wronskian of 
f q and its complex conjugate. However, the asymptotic regime rj 3> 1, and x fixed is 
completely determined by the scaling form of the solutions given by eqs. ( |3.14| ), ( |3.15| ) 



with the modulus squared of the coefficient A determined by the sum rule and given 
by 



In the asymptotic regime rj ^> 1, x fixed which is dominated by the scaling solution 
we see that only the long wavelength modes are relevant, i.e. q = x/t] <^ 1 (in units 
of mfj), justifying keeping only the classical part J c i in eq. (|2.44|) and using the cutoff 
k ~ rriR in the integral in Q2.45| ). 

The amplitude of the long wavelength modes q = x/rj that dominate the scaling regime 
become non-perturbatively large as can be seen in figs. d^,§) and is a consequence of the 
sum rule ( |3.8| ) that fixes the amplitude of the coefficient of the scaling solution. The 
fact that the long-wavelength modes become non-perturbatively large and dominate 
the dynamics is a consequence of the early time linear instabilities that result in an 
exponential growth of these modes. 

A noteworthy feature of the scaling solution is that its phase has frozen i.e. it became 
time independent and completely determined by the phase of Aq. We note that if 
the two linearly independent solutions in eq. (|3.12|j3.13l ) had the same amplitude the 
phase of the mode functions will be time and momentum dependent. The freezing 
of the phase and the non-perturbative large amplitude of the long-wavelength modes 
entail that these mode functions that originally had quantum initial conditions (as can 
be seen by restoring h to obtain the quantum commutators) had become classical. 



B. Scaling corrections 



Although the leading asymptotic behavior of the mode functions is determined by the 
scaling forms (|3.14j , [3.15| ) the subleading contributions determined by the terms containing 
the Neumann functions N p+ i(x) in eqs. ( |3.12 , 3.12| ) are important corrections to scaling 
behavior and are required for the self-consistency of the solutions. 

The analytic and numerical result C 2 (r])M 2 (rj) v — — ( ' 2p+ '^ 2p+1 ^ entails two very strin- 
gent constraints. The first one leads to the cancellation between the tree level contribution 
and the backreaction term g£(r?) = 1 to leading order in C(r]) and leads to the sum rule eq. 
(ET5|). However the fact that C 2 (ri)M 2 ( V ) ?? = 1 - (2p+ g 2p+1) 
T e ff(r])/T c <C 1 implies that 



in the asymptotic region when 



C 2 ( V )M 2 (rj)=g q 2 alq | f q (r)) \ 2 - C 2 {r]) 



??>i 



(2p + 3)(2p 
4?7 2 



(3.17) 



Inserting the solutions ( |3.12|) - (|3~T3|) in eq. (|3.17|) , changing integration variables to q = x/rj 
and taking the asymptotic limit rj 3> 1, x fixed inside the integral it is seen that the leading 
contribution arises from the terms proportional to J p +\(x) which cancel C 2 {r]). The next 
order term arises from the crossed term between the contribution proportional to the J p+ i (x) 
and that from the N p+ i(x) which leads to another sum rule 
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2sWa ix JM*) W-) - - <2P + ? ( f +1) (3.18) 



O / ^Ti\~/ i'Ti \ ~ / ^1 O 

77^ jo 4 r\ z 

which describes the asymptotic behavior of C 2 (rj)M 2 (rj) 



Using the result J24 



oc 



X 



dx J p+ x(x) N p+ i(x) 



(2p + 3)(2p+l) 



10 lb 
we find that the sum-rule ( |3.18| ) constrains the value of Re[v4 -Bo] as follows, 

Re[A B*\ = -- 

thus we are led to the conclusion that the self-consistency condition indeed requires a non- 
vanishing correction to scaling. 

Eq. ( |3.18[) reveals that B oc 1/y/g therefore of the same order as Aq. Thus we conclude 



that although the asymptotic behavior of the mode functions is dominated by the scaling 
form, the corrections to scaling embodied in the subleading contributions are very important 
for the self-consistency of the solution, they have non-perturbative amplitudes of Q (1/y/g) 
and as argued above the coefficients B q do not vanish because of the Wronskian condition 
on the mode functions. 



C. Classicality 

We argued above that the long-wavelength modes become classical in the sense that 
their amplitudes become non-perturbatively large 0(1/ y/g) and their phases freeze out. An 
important bonus of studying the dynamics in terms of a quantum density matrix is that the 
classicalization of long-wavelength fluctuations can be quantified in terms of the probability 
density (functional) in field space. Just as in quantum mechanics the probability density is 
the diagonal density matrix element in the Schrodinger representation, i.e. 

V($) = p[$ , $, t] = ]J A4(t) exp {-[A k>R (t) + B k (t)} $ fc ■ $_ fc } (3.19) 
k 

uik/Ti and in terms 



For long- wavelength modes k << T{ we can approximate 2tanh 
of conformal time and the mode functions that obey the equations fl2.46|) -( f2.50|) we find 



m ^n^)exp{-q^g^} ,3,0, 

Thus at any given fixed (conformal) time, configurations with amplitude $^ ~ 
VT~i\fk(v)\/C(r]) (in units of tur) have a probability ~ 0(1) of being represented in the 
statistical ensemble. Since Tj > T c w 1/g and in the asymptotic regime (7] 3> 1 ; x fixed) 
\fk(v)\ ^-/\/9i these are large amplitude 0(1/ g), long wavelength k < l/rj field configura- 
tions. 
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Thus whereas asymptotically for r] ^> rjj^L these large amplitud long- wavelength con- 
figurations are represented in the ensemble with probability 0(1), in the initial density 

matrix with |/fc(0)| ~ 1 these field configurations are represented in the ensemble with 

_ i 

probability a e » <C 1 in the weak coupling limit. Under time evolution the Gaussian 
probability density (functional) for long-wavelength modes spreads out and large amplitude 
long-wavelength configurations acquire non-vanishing probabilities of being represented in 
the statistical ensemble. Thus in the asymptotic scaling regime, a typical field configuration 
found in the statistical ensemble will have amplitude 0(1/ g) and the Fourier transform of 
its spatial profile will be dominated by long- wavelength modes k <C 1/r]. This observation 
leads to a semiclassical stochastic description in terms of semiclassical field configurations 
that describe a typical member of the ensemble for t] ^> t)nl 

^)~E cos ^ • f + $ ( 3 - 21 ) 

where the phases 52 are stochastic with a Gaussian distribution in order to reproduce the 
field correlation function obtained from the Gaussian quantum density matrix given by the 



integrand in ( 2.44Q . The phases 5k represent the phases of the coefficients A q in the scaling 



solutions ( |3.14| )-( |3~T5|) which as argued above is time independent. We note that a spatial 
translation can be absorbed into a redefinition of the phase 5k and the stochastic nature of 
this variable in terms of a Gaussian probability distribution restores translational invariance 
in the ensemble averages of the semiclassical but stochastic field configurations. 



IV. CONSEQUENCES OF SCALING: DYNAMICAL CORRELATION LENGTH 

AND EQUATION OF STATE 

A. Dynamical correlation length 

As described in the introduction in known systems the dynamics of phase ordering leads 
to the emergence of a dynamical correlation length £(t) which determines the size of the 
correlated regions ||- 0. This dynamical correlation length plays the same role in the 
description of the dynamics as the static correlation length does in static critical phenomena. 
When the correlation length becomes much larger than the typical microscopic length scale 
in the system it becomes the only relevant length scale and much in the same way as in 
static critical phenomena a scaling regime emerges where the correlation length provides the 
natural scale for all dimensionful quantities. The emergence of the dynamical correlation 
length in dynamical critical phenomena is revealed by the equal time correlation function 
($ a (s, i)$ 6 (0, t)). From the discussion in the previous sections, this correlation function 
will be dominated by the long-wavelength modes, hence we use the high temperature limit 
Tj >> k to find 

(*»( 2 , t «o, i )>^-r,/^!^f^! (4,, 

In the scaling regime when the mode functions are of the form (|3.6|) and multiplying by the 
coupling constant A to write the result in a more familiar manner we find 
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\(<S> a (x,t)<S>\d,t))=5 a ' b gD(z) 

\ r°° t 
D(z) = — xdx sin2xz \F{x)\ 2 ; z = — (4.2) 
2z Jo 2r\ 

where we have introduced the scaling ratio z = r/2r]. Figures ||, |l^ show gzD( 
function of z for rj = 240, 400 (in units of tor) respectively for radiation dominated cosmology 
obtained from the integration of the mode functions fk{v)> an d figure 11 shows gzD(z) for 
the matter dominated case for r\ = 400. A remarkable feature of the correlation function 
D(z) is that it becomes of the order g for r > 2i] (i.e, z > 1) in both cases, obviously as a 
result of causality. This result leads to the conclusion that the dynamics of phase ordering 
is described in the scaling regime as the growth of correlated regions of comoving size 2rj, 
the reason for the factor 2 is that one edge of this region is localized at the point and the 
other at r and the boundary of this correlated region recedes at the speed of light. 

Hence we obtain that the comoving dynamical length scale that determines the typical 
(comoving) size of a correlated domain is given by 

icomijl) = V (4-3) 

The physical dynamical correlation length is given by 

W(^) = C(V) V = d H (v) ( 4 - 4 ) 

with dii(v) the size of the causal horizon. Hence we conclude that the dynamical correlation 
length is exactly the size of the causal horizon and within one horizon there is exactly one 
correlated region within which ($ • $) « m 2 R /\n i.e, the mean square root fluctuation of 
the field is probing the broken symmetry ground state. This result has been obtained via 
the full numerical evolution of the mode functions but can be understood analytically as a 
consequence of the scaling solutions ( |3.14| , |3.15| ) that dominate the asymptotic regime for 
long-wavelengths. Replacing J-{x) in (|4.2|) by the scaling form of the solutions given by 
Q3.14| )-( p.l5f) in the radiation and matter dominated cases we find 



z") J F\h, I; 4; 1 — 2; 2 16(1 — z) for Radiation dominated 



v 2\3 z?ri 5. 

1.2' 2' 

gzD(z) = { (4.5) 
(1 - z 2 f F[\, f ; 6; 1 - z 2 }@(l - z) for Matter dominated 

where F[a, b; c; z] is the hypergeometric function. A numerical evaluation of these expressions 
agrees to a very high level of precision with the result obtained above from the numerical 
integration. 

Moreover, these particular hypergeometric functions can be expressed in terms of com- 
plete elliptic integrals |25]| 

15 o-, 32 

-, -;4; 1 - z 2 } = 

L 2'2' ' J 15tt 



1 - z 2 f F[-, -; 4; 1 - z 2 ] = — z 2 {z 2 - 9)K(VT^) + (3 + 7z 2 - 2z 4 )E(VT^ 2 ) 



■ 17.. „ 512 



[l - z 2 f F[-, ^ 6; 1 - z 2 } = (1 - z 2 ) 2 [(45/ - 7z 6 + 22z A - llz 2 + 15)ir(v / I^) 



;i5z 8 + 31z 6 + 12/ - llz 2 + 23)£(v / T^2) , (< 
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where K(k) and E(k) are complete elliptic integrals of the first and second kind, respectively. 

The power spectrum in the asymptotic scaling regime is dominated by the long- 
wavelength modes and is given by 



9 



(hi) 



p+i 



(4.7) 



where p — 1, 2 for radiation dominated or matter dominated, respectively and where Z is a 
constant of order one. At long times rj ^> 1 this power spectrum becomes strongly peaked 
at k = and receives contribution only from a narrow region k < 1/rj as evidenced by figs. 
(§-0). Since the sum rule constrains the total integral of the power spectrum to be 



je dkS(k,t) = constant 



(4.8) 



we are led to the conclusion that at asymptotically long times the power spectrum is sharply 
peaked at k = becoming asymptotically a delta function. This is the signal of the formation 
of a zero momentum condensate, much in the same manner as observed in the case of 
Minkowski space-time P,[i~T|| and also in dynamical critical phenomena in condensed matter 
systems [|]|§. 



B. Equation of state 

The expectation value of the energy momentum tensor in the non-equilibrium density 
matrix is of the fluid form in terms of the energy density e and pressure P given by 



A - I 

N £ ~ 2 



o l 2N 1 



) 


+ - [ d 3 k coth 








2.1 







\Mt)\ 2 + w2(t)\Mt)f 



A(p + £ ) = \Jd 3 k coth 







l2T t _ 





3a 2 (t) 



(4.9) 
(4.10) 



where (pk(t) obeys the equations of motion ( |2.24| ) with W%(t) given by eq. ( |2.13j ) and we 
had set to zero the coupling to the Ricci scalar. Using the equations of motion Q2.24 ) and 
the definition of the self-consistent frequencies W%(t) ( |2.13| ), it can be easily verified that e 
and P satisfy the covariant conservation law 



e(t) + 3 - [e(t) + P(t)} = 



(4.11) 



Passing on to conformally rescaled mode functions and conformal time and writing the 
asymptotic forms of the scale factor in radiation dominated and matter dominated domi- 
nated cosmologies as 



with p = 1 for Radiation dominated ; p = 2 for Matter dominated (4.12) 



we will focus our discussion on the behavior of the pressure and the energy density in the 
scaling regime wherein the conformally rescaled mode functions in conformal time are given 
by eq. ( |3.14| , p,15|) , which can be handily written as 



20 



h(t) -AoV xP+1 



(4.13) 



in either case. The sum rule ( |3.10|) with vanishing coupling to the Ricci scalar leads to the 
following identity in the asymptotic regime 



c 2 (v) 



2 A ,£ ^ C"(V) 



(2p + 3)(2p 
A-rf 



+ 0(l/r] 4 



(4.14) 



up to perturbatively small corrections. The leading term in this sum rule is given by the 
sum rule ( |3~5|) which upon using the asymptotic scaling mode functions ( p.l4j )-( |3T5f) can be 
written in the following compact form 



\A \ 2 g J dx 



xp 



-2p 



(4.15) 



Using this result and after some straightforward algebra we find that the asymptotic behavior 
of the energy density and pressure are given by 



A 



\A \ 2 g2 2 P 
C 2 (r]) rj 2 



x dx< 



1 J p+1 (x) J' p+1 (x) 
^2~ v > — + 



1 2 



xf> 



+ x z 



J, 



p+1 



X 



1 2 



XP^ 



(4.16) 



N 



e + P) 



\A \ 2 g2^ 
C 2 {rj) rf 



x dx< 



1 2 



X 



■p+1 



XP 



+ 



X 



1 2 



X- 



•p+1 



(4.17) 



where the prime stands for derivative with respect to the argument. The integral is inde- 
pendent of conformal time and the only dependence on r\ is in the prefactors, the product 



C 2 ( V ) rf oc 



thus we find that in either case 



C A {rj) for Radiation dominated 
C 3 (rj) for Matter dominated 

£ OC Pback{rj) 



(4.18) 



(4.19) 



where pback(v) is the background energy density in either radiation dominated (pback(v) (X 
C~ A (r])) or matter dominated (pbackiv) ^ C ^ 7 ?))- Since e and P satisfy the covariant 
conservation equation and so do the background energy density and pressure we conclude 
that the fluid resulting from the fluctuations of the scalar field in the asymptotic scaling 
regime obeys the same equation of state as the background fluid. Indeed careful evaluation 
of the integrals using the properties of the Gamma function ]24| as analytic functions of the 
variable p lead to the following remarkable results 



A 

N 



e{ri) C 2 {n) rj 2 = 3 p 



p + 1/4 
p — 1 



(4.20) 
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±[e(r,) + P( V )] C 2 ( V ) = 2 (p + 1) (4.21) 
A/ p — 1 

leading to one of the important results of this article 

Q -, ' -) \ i | for Radiation dominated (p = 1) 
' (4.22) 

for Matter dominated (p = 2) 




The simple poles at p = 1 in eqs. (|4.20| )- (|4.21| ) reflect ultraviolet logarithmic singularities 



in the integrals which have been evaluated with an upper limit taken to infinity. However 
we remark that the scaling form requires that the upper limit be of order x <C 1 /rj so that 
the coefficient of the Bessel functions can be replaced by momentum independent constants 
leading to true scaling solutions. A numerical evaluation of the integrand in ( |4.16| )-( fLl7| ) 



reveals that these are strongly peaked near x = and most of the contribution arises from 
the interval x < 5 even for p = 1. The results ( 4.2U| )- (|1.21j) are the analytic continuation 



of these integrals in the variable p taking the upper cutoff to infinity. The variable p thus 
acts as a regulator much in the same way as analytic or dimensional regularization and the 
equation of state is insensitive to the regularization procedure. 

Defining e = 5p with e the energy density given by eq. ( |4.16| ) and 5p as the contribution to 
the energy density from the fluctuations of the scalar field, the fact that Sp/pi, ac k =constant 
i.e, independent of time is the statement that the power spectrum for the density fluctuations 
is of the Harrison- Zeldovich form, i.e., scale invariant, for long-wavelength scalar density 
perturbations [l].§]. 



C. Universality of scaling 

Although we have studied radiation and matter dominated cosmologies separately while 
the most physical scenario involves a smooth transition between the two regimes we now 
argue that if the transition between the two regimes occurs over a very long time scale our 
analysis leads to the conclusion that scaling will be a robust feature in the regimes in which 
the scale factor is dominated by a power law. In the asymptotic regime, when the effective 
time dependent mass vanishes, the sum rule 

9 f Q 2 dq \ f q (rj) \ 2 =C 2 ( V ) (4.23) 
Jo 

is fulfilled for an arbitrary scale factor, not necessarily a power law, however if C(rj) is not 
a power law, the solution is not of the scaling form. 

Consider the scale factor in the general case in which the background fluid has a radiation 
and a matter component as given by ( |2.52| ),i.e, 

C( V ) = 1+H irj + Jf-^- V 2 (4.24) 

4(1 + n) 

in a realistic scenario the ratio of the initial matter density to the initial radiation density 
is Ti <C 1. Therefore in the regime — ^> Hi rj 3> 1 the scale factor is approximately that 
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of a radiation dominated cosmology and therefore a power law, hence the sum rule above 
( |4.23p leads to the scaling form of the solution. During the time scales of the crossover 
between radiation and matter domination r^if,?] « 1 the scale factor is not a pure power 
law. The sum rule above still holds but it does not entail a scaling form for the mode 
functions. However for H^Ti 3> 1 again the scale factor is a power law now describing a 
matter dominated cosmology and the sum-rule again leads to a scaling form for the solution 
in the long- wavelength limit. However when the new scaling form emerges the solution does 
not depend on the past history, i.e. on the violations of scaling during the crossover or the 
previous scaling solution in the radiation dominated regime. The reason for this is that 
the sum rule above is local in time and fixes the absolute value of the coefficient Aq and 
the scaling form in terms of the Bessel function, however the phase of the coefficients A q 
and certainly of Aq will depend on the history and will not be determined by the sum-rule. 
However in the semiclassical stochastic description advocated in the previous section this 
phase is treated as a stochastic variable with a Gaussian distribution function. Hence all of 
the initial information and the history contained in the phases of the scaling solution can 
be treated stochastically for the long- wavelength components. 



V. COMPARISON WITH THE O(N) NON LINEAR SIGMA MODEL: 

There are many similarities but also many differences with the results obtained in the 



case of the classical non-linear sigma model in refs. flgjlofl . The non-linear sigma model 
describes the non-linear interactions of Goldstone bosons in a broken symmetry phase. To 
begin with let us establish an important difference: the non-linear sigma model cannot 
be used to describe the phase transition because the fields are constrained to the vacuum 
manifold. Hence all of the details of the dynamics of the phase transition, the linear in- 
stabilities, growth of long-wavelength modes, the different time scales and regimes and the 
explanation of why the power spectrum is peaked at long-wavelength as a consequence of 
the early exponential growth of long-wavelength fluctuations simply cannot be captured by 
the non-linear sigma model. These details that are deeply dependent on the details of the 
phase transition have important consequences: the amplitude of the mode functions in the 
scaling regime is non-perturbatively large in the weak coupling limit since Aq oc 1/ J~g. This 
amplitude is not arbitrary, if h is restored, the initial amplitudes of the mode functions are 
ofO{V%), the non-perturbative growth of the amplitude of long-wavelength fluctuations is 
a consequence of the linear exponential instabilities which are a dynamical hallmark of the 
phase transition. Furthermore the corrections to scaling embodied in the contribution from 



the Neumann functions in eqs.( |3.12"D -( |3.13D have very precise coefficients B q since these are 



constrained by the Wronskian condition and the self- consistency condition ( 3.18|) . Both the 
coefficients A q ; B q carry information from the initial conditions because of the constancy 
of the Wronskian, the long wavelength limit Aq ; Bq are of order 0(1/ ^/g) as a consequence 
of the self-consistency conditions fl3.16|) -( [3.18D and the phases can only be determined from 



the numerical evolution. 

There are also very important similarities: in the asymptotic regime the scaling form of 
the solution ( |3.14| )-( |3.15| ) is the same in the model studied here and in the O(N) non-linear 
sigma model |15|JT6[1 . In the non-linear sigma model the coefficient of the Bessel function is 
fixed by the length of the O(N) vector, this length is interpreted as the vacuum expectation 
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value of the field. In fact in the model studied here the sum-rule Q3.5Q can be identified with 
the same constraint. The similarity of the asymptotic scaling solutions is a consequence of 
the fact that eq.(|3.5|) is local as discussed above. 

The description of the non-equilibrium dynamics via the evolution of the quantum den- 
sity matrix leads to a clear interpretation of the emergence of a semiclassical stochastic 
description of long-wavelength fluctuations at long times. Our analysis thus provides a 
consistent microscopic, quantum field theoretical derivation for the classical stochastic treat- 



ment of the Goldstone modes in the asymptotic region which was used in references [|15l , |16 
Furthermore, our detailed analysis quantifies the time scales for which such semiclassical 
stochastic description is valid and goes beyond providing corrections to scaling. 

In ref. fl5j the energy density of the scalar field in the asymptotic scaling regime was found 
to be proportional to that of the background in the matter dominated case and differing from 
that of the background by a weak logarithmic dependence on r] in the radiation dominated 
regime. The logarithmic dependence has a counterpart in our result in the form of the 
pole in the energy density and pressure as p — > 1 in the expressions ( |4.20| )-( fOl] ) obtained 
via an analytic continuation in the variable p, as we pointed out the equation of state is 
independent of the regularization. Furthermore we emphasized that both the energy density 
and the pressure are dominated by modes k << l/rj i.e. superhorizon modes. 



VI. CONCLUSIONS AND FURTHER QUESTIONS 

In this article we have studied in detail the non-equilibrium dynamics of a symmetry 
breaking phase transition in a spatially flat radiation and or matter dominated FRW back- 
ground. Anticipating the necessity for a non-perturbative treatment, we studied the linear 
sigma model with a scalar field in the vector representation of O(N) in leading order in 
the large N limit. The Liouville equation for the quantum density matrix is solved in this 
limit from which we extract the necessary correlation functions to leading order. The advan- 
tage of working with a density matrix is that this description allows a clear interpretation 
of the emergence of a semiclassical description in terms of a field probability distribution 
functional. 

The main goal of our study is to provide a thorough understanding the process of phase 
ordering beginning from a state of local thermodynamic equilibrium at an initial tempera- 
ture larger than the critical, the cosmological expansion triggers the phase transition when 
the effective time dependent temperature falls below the critical. The key issues that we 
address in this article are the emergence of a scaling regime and of a characteristic dynami- 
cal correlation length that determines the spatial extent of the correlated regions of broken 
symmetry and the consequences of such scaling behavior upon the equation of state of the 
field fluctuations. 

The non-equilibrium dynamics after the phase transition is characterized by three dis- 
tinct time scales: during the early stage after the phase transition the dynamics is dominated 
by the exponential growth of long- wavelength fluctuations associated with spinodal insta- 
bilities. This is essentially a linear regime in which the backreaction can be ignored for 
weak self-coupling and an analytic description is available. An intermediate time scale is 
defined when the backreaction from the self-consistent mean-field is comparable to the tree 
level terms in the equations of motion and determines the onset of non-perturbative and 
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non-linear dynamics which must be studied numerically. A third, asymptotic time scale 
reveals the emergence of a scaling regime for radiation or matter dominated FRW back- 
grounds. This stage is dominated by large amplitude fluctuations with wavelengths of order 
of or larger than the causal horizon. This regime is characterized by a dynamical physical 
correlation length £ p h ys = dn{t) with du{t) the size of the causal horizon and the onset of 
a non-equilibrium condensate at zero momentum. The dynamical correlation length deter- 
mines the size of the correlated domains inside which the field fluctuations probe the broken 
symmetry states, hence there is exactly one correlated domain per causal horizon. The field 
correlations vanish for distances larger than djj{t) by causality. In this regime the phases of 
the long-wavelength quantum mode functions become time independent and their amplitude 
becomes non-perturbatively large, the approach via the density matrix reveals the emergence 
of a semiclassical but stochastic description and provides a microscopic justification for a 
semiclassical stochastic treatment in the asymptotic regime. 

A remarkable corollary of the scaling solution of the equations of motion is that in the 
asymptotic regime the equation of state of the fluid described by the fluctuations of the scalar 
field is the same as that of the background fluid that drives the dynamics of the scale factor. 
An important consequence of this behavior of the fluid is that Sp/pb ac k =constant which 
results in a Harrison-Zeldovich spectrum of scalar density perturbations for superhorizon 
wavelengths 0)0- 

The self-consistency of the equations of motion entails very precise corrections to scaling. 
We argue that the scaling solution is a universal feature of a scale factor that is a power law 
and is independent of the crossover between the radiation and matter dominated regimes. 

Some important aspects remain to be explored further. A thorough investigation of 



cosmological perturbations has been performed in ref. [fLq| . An important aspect of the 
models under consideration is that the perturbations in the energy momentum tensor are 
non-linear in terms of the field fluctuations, unlike the case of scalar field perturbations 
in the inflationary stage. The correction to the energy density is given by the expectation 
value of the energy momentum tensor of the fluctuations and even to leading order in the 
large N limit, this is a quadratic form in terms of the fluctuations [see eq. (|4.9| )1 . This 
feature has an important consequence in the spectrum of primordial density perturbations: 
decoherence effects tend to suppress the acoustic peaks at large angular momentum I (small 
angular resolution) (see JTBJ for a clear treatment in the non-linear sigma model). As we 



have argued above, there are very precise and important corrections to scaling that are 
manifest for subhorizon wavelengths, i.e. k ^> 1/n. An important possibility is that these 
corrections to scaling affect the correlations of the energy momentum tensor and therefore 
the power spectrum of scalar density perturbations on small angular scales. This would 
arise from the interference effects between the scaling contributions and those of the scaling 
corrections, much in the same manner as the next-to-leading contribution in the sum rules 
discussed above. As we have highlighted these interference effects are necessary for the 
self-consistency of the method leading to a precise form of the scaling corrections. These 
corrections had not been taken into account in the analysis of the power spectrum at small 
angular scales and whether they lead to significant corrections and or a dramatic change of 
the picture of decoherence effects in causal perturbations is an open question that merits 
careful study. 

Another important aspect that requires further understanding is the description of the 
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non-equilibrium dynamics in the case of a single scalar field which cannot be addressed 
reliably in the large N expansion. In the large N limit an important feature of the phase 
diagram of the theory is that the coexistence and the spinodal line merge as a consequence of 
the Ward identity associated with the continuous symmetry |l3|Jl^| and the long-wavelength 
physics described by the dynamics of Goldstone bosons. Static renormalization group ar- 
guments in scalar field theories lead to a Maxwell constructed free energy with a region of 
phase coexistence that joins the minima of the free energy and no metastable region 
for all values of N with no basic difference between the cases of a single scalar field or a mul- 
tiplet of fields. It is important to find a consistent non-perturbative approximation scheme 
for iV = 1 to study the dynamics in this case. Another important aspect to be studied 
further is whether the scaling solution survives 1/N corrections, in particular it is clear that 
in next to leading order in the large N approximation another time scale associated with 
collisional processes should emerge and the relevant question is whether this microscopic 
scale will modify substantially the scaling solution and in particular the equation of state. 
These and other related questions are currently under consideration. 
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FIG. 8. (-4/15)rj 2 C 2 (rj)M 2 (ri) vs. r\ (in units of m^ 1 for J = 3, g = 1(T 5 . Radiation 
dominated universe. 
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